Binary self-similar one-dimensional quasilattices: 
Mutual local-derivability classification and substitution rules 



Masashi Torikai 
Department of Physics, Kyushu University, Fukuoka 812-8581 

Komajiro Niizeki 

Department of Physics, Graduate School of Science, Tohoku University, Sendai 980-8578 

Takashi Odagaki 
Department of Physics, Kyushu University, Fukuoka 812-8581 

Self-similar binary one-dimensional (ID) quasilattices (QLs) are classified into mutual local- 
' derivability (MLD) classes. It is shown that the MLD classification is closely related to the number- 

theoretical classification of parameters which specify the self-similar binary ID QLs. An algorithm 
£NJ ■ to derive an explicit substitution rule, which prescribes the transformation of a QL into another QL 

in the same MLD class, is presented. An explicit inflation rule, which prescribes the transformation 
of the self-similar ID QL into itself, is obtained as a composition of the explicit substitution rules. 
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1 !• INTRODUCTION 

(N 

A quasiperiodic structure lacks the periodicity but has a long-range positional order. Quasiperiodic structures have 
been intensively investigated since the experimental discovery of quasicrystals, which have quasiperiodic structures 
with non-crystallographic rotational symmetries!] In order to discuss the peculiar properties of quasicrystals, many 
models for the structures of quasicrystals have been proposed. The most successful model is based on a quasilattice 
^ ■ (QL), which is a quasiperiodic and discrete set of points. 

The mutual local-derivability (MLD) is one of the most important classification schemes for QLs. Em According to a 
recent study of a tight-binding model on a one-dimensional (ID) QL, two ID QLs belonging to different MLD classes 
belong to different universality classes with respect to their one-electron properties .□ 

Some of the quasiperiodic structures have self-similarities. Jn particular, all the quasicrystals found experimentally 
are closely related to self-similar quasiperiodic structures .tU'El On the other hand, the self-similarity has been used 
£j , successfully for a rigorous study, of the one-electron properties of ID QLs because it matches well renormalization 
J* ■ methods such as the trace map.tTtl 
. £h ! Many studies have focused on self-similar ID quasiperiodic structures, but some important points have been left. 
Kaneko and Odagaki have discussed the quasiperiodic sequence of letters (i.e., objects without geometricaL|k|ngth), 
which are produced by associating different letters to different types of intervals (atomic distance) in QLs.MLil They 
have shown a procedure to obtain an explicit inflation rule (IR) , which prescribes the transformation of a self-similar 
quasiperiodic structure into itself, for a binary quasiperiodic sequence of letters characterized by a quadratic irrational. 
However, since the quasiperiodic sequence of letters is not a geometrical object, the geometrical meaning of the explicit 
IR is obscure. Binary sequences derived from substitution rules (IRs in our terminology) are extensively investigated 
in Ref. |^. They found that there exist hierarchical series of substitution rules which have good properties with 
respect to the trace map. They discussed interrelations among different substitution rules. This work also treats the 
binary sequences as sequences of letters. Yamashita has investigated the relation between the projection method and 
two series of binary quasiperiodic sequences of letters with inflation symmetries.E3 However, he has not treated the 
sequences as geometrical objects. As a consequence of considering the quasiperiodic sequence of letters instead of the 
QL itself, these works have not revealed the geometrical meaning of the self-similarity and the IR . 

In order to consider the geometrical meaning of the self-similarity and the IR, it is necessary to investigate the QLs 
instead of the sequences of letters. In Ref. |l3|, the self-similarity of ID QLs has been investigated, and the geometrical 
meaning of the self-similarity has been revealed. However, a general algorithm to obtain an explicit IR has not been 
shown. In Ref. [l4], the inflation symmetry of binary and ternary ID QLs has been extensively investigated but this 
study has been restricted to the cases where the quadratic irrationals characterizing the QLs belong to two infinite 
series, one of which is the series composed of all the quadratic irrationals of rank one in our terminology to be presented 
in a later section, while the other is a subseries of quadratic irrationals of rank two. 
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In the present paper, we restrict our considerations to self-similar binary ID QLs but we do not impose any condition 
on the quadratic irrationals characterizing the QLs. The purpose of the present paper is to classify all the self-similar 
binary ID QLs into MLD classes and to derive explicit substitution rules (SR), which relate different members of a 
single MLD class. 

In the next section, we make a review of the self-similarity of binary QLs. Here we identify parameters which 



characterize these QLs. In Sec. we classify the QLs into MLD classes on the basis of the number-theoretical 
properties of the parameters. Also, we give an explicit SR which combines two QLs in a single MLD class. In Sec. 
IV, we discuss composite SRs and IRs. We show several examples of self-similar binary Q Ls in Sec. |v|. In Sec. VI 



we discuss symmetric SRs and symmetric IRs. Miscellaneous subjects are discussed in Sec. VII. In the last section, 
we summarize the results, and discuss a few remained subjects. 

In order to develop our theory, several results and technical terms in the number theory of quadratic irrationals 
are necessary, and they are explained in Appendix |A|. We will sometimes use them in later sections without explicit 
references to the appendix. We suggest the readers jumping to Appendix at this point although it is considerably 
long. 

Finally, a list of some of abbreviated words appearing in the present paper is presented: inflation rule (IR), 
substitution rule (SR), mutual local- derivability (MLD), unimodular Frobenius matrix (UFM), modular equivalent class 
(MEC), quasi-permutation-matrix (QPM), reduced quadratic irrational (RQI), and quasi-reduced quadratic irrational 
(QRQI). The last five words are defined and discussed in Appendix. 



II. SELF-SIMILARITY OF BINARY ID QLS 
A. Binary QLs in the projection method 

A QL is derived, as shown in Fig. [[], with the projection method from a 2D periodic lattice A, which is called a 
mother lattice. The 2D space embedding A is spanned by two orthogonal subspaces E\\ and E±, which are called 
the physical space and the internal space, respectively. We shall use throughout this paper the coordinate system in 
which the x axis and the y axis are chosen to be E\\ and E±, respectively. The physical space is assumed to be in 
an incommensurate configuration with A. One of the two basis vectors, a and b, of A is assumed to be in the first 
quadrant and the other in the fourth quadrant, so that the fundamental parallelogram which is spanned by a and b 
is cut by . If all the lattice points within a strip which is parallel to E\\ are projected onto Eu , we obtain a QL as 
shown in Fig. [j]. The QL divides naturally Eu into intervals, so that we shall identify the QL with the relevant ID 
tiling of Eh. In the present paper, we restrict our considerations to binary QLs. 

We call the section of the strip cut by E± a window. A binary QL is obtained as shown in Fig. [j] if the size of the 
window is equal to the vertical width of the fundamental parallelogram. The two types of intervals are denoted as a 
and f3, whose lengths, a and 6, are equal to the x-components of the two basis vectors. It is essential in obtaining 
a binary QL that one of the two basis vectors of A is in the first quadrant and the other in the fourth quadrant 
and that the size of the window is equal to the vertical width of the fundamental parallelogram. Therefore, we shall 
call such a set of basis vectors a canonical set and the relevant window a canonical window; a canonical window is 
associated with a canonical set. Since there exist an infinite variety of canonical sets of basis vectors, there exist an 
infinite variety of canonical windows. This does not mean, however, that there exist an infinite variety of binary QLs 
with a common mother lattice because different QLs can be (geometrically) similar to one another. It will be shown 
in a later section that, under a certain condition, there exist only a finite number of binary QLs which are derived 
from a common mother lattice but not mutually similar; the condition is that the mother lattice has the hyper scaling 
symmetryX.3 

Two QLs are called to be mutually locally isomorphic (LI) if every finite part of one of them is always found in 
the other and vice versa; two QLs being mutually LI cannot be distinguished macroscopically from one another. A 
necessary and sufficient condition for two QLs to be mutually LI is that the relevant two windows have a common 
size. In the present paper, we will not distinguish two QLs if they are mutually LI. We denote Q ~ Q' when two QLs, 
Q and Q', are mutually LI. A QL derived from a mother lattice is uniquely determined by the window, W, and the 
QL specified by W is denoted as Q(W). We will not hereafter distinguish the window, W, and its size, \W\. 

It is known that a Fourier module, i.e., an additive group formed of reciprocal lattice vectors of a QL, is given 
by projecting the Fourier module of the mother lattice onto the reciprocal physical space. Therefore different QLs 
derived from a single mother lattice but with windows of different sizes have a common Fourier module. 

The scale of the internal space, E±, is irrelevant on the construction of a QL because the projection is made along 
E± in the projection method. Hence, two mother lattices which only differ in the scale of the internal space are 
considered to be isomorphic, and are not mutually distinguished. On the other hand, a scale change of the physical 
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space, E\\, results in a scaling of a QL, so that to fix the scale is equivalent to fix the lattice constant of the QL. In 
a usual classification scheme of QLs, two QLs being geometrically similar to each other are assumed to be identical. 
Then, two mother lattices which differ in the scales of the physical space and/or the internal space are sometimes 
assumed to be identical. 



B. Self-similar binary QLs 



The Fibonacci lattice is a representative self-similar binary QL. The sizes of two types of intervals, a and /3, of the 
Fibonacci lattice are equal to 1 and tq, respectively, where tg is the golden ratio: tq := (1 + y/o)/2. If we substitute 
the intervals /3 and a(3 for a and (3 in the Fibonacci lattice, respectively, we obtain again the Fibonacci lattice. This 
transformation of the Fibonacci lattice into itself is known as the inflation rule (IR), which is represented as 



(3, /3-> ap. 



(2.1) 



Each interval is scaled up by tq on the inflation. 

The self-similarity of a QL is represented generally by an IR which prescribes how the inflated versions of the two 
intervals are divided into the original intervals £3 The IR acc omp anies naturally the inflation matrix M, which is a 
unimodular Frobenius matrix (UFM) as defined in Appendix A 2. For example, the inflation matrix of the Fibonacci 
lattice is given by 



M 



1 

1 1 



(2.2) 



The golden ratio, tg, is not only the Frobenius eigenvalue of M but also the companion quadratic irrational of M. 
Since the explicit form of the inflation matrix M of a self-similar binary QL depends on the order of the two intervals 
in the defining equation for M, we shall fix the ord er so that the first interval is the shorter of the two. Then, M 

A self-similarity of a QL is derived in the projection method 
this symmetry is an automorphism of A and leaves the 



is a quasi-normal UFM as defined in Appendix A 2 



only when the mother lattice has the hyperscalingjSymmetry: 
physical space and the internal space invariant EjEj 

Since the inflation matrix M is unimodular, it determines a 2D lattice A with the hyperscaling symmetry as shown 
in Appendix A 1 . The physical space (or the internal space) is chosen to coinci de w ith the expansive (or contractive) 
principal axis of the linear transformation T associated with M (see Appendix A 1 ) . More precisely, T enlarges along 
E\\ by r, i.e., the Frobenius eigenvalue of M, but shrinks along E± by r _1 = |f|. The left Frobenius eigenvector (a b) 
with a, b > satisfies the equation 

r(o b) = (a b)M, 

while the second left eigenvector, (a± b±), satisfies the equation 

f(a ± b±) = (a± b±)M. 

The two basis vectors of A are given by 

' b 



b = 



(2.3) 
(2.4) 

(2.5) 



It follows from Eqs. ( p.3j ) and (2.4) that T is an automorphism of A; T yields the hyperscaling symmetry of A. Since 
the row vector (a± b±) is not a Frobenius eigenvector, a±_ and b± must have different signs. We can assume that 
a± < and bj_ > 0. Then, a and b are in the fourth quadrant and the first quadrant, respectively, so that {a, b} is 
a canonical set. The QL Q(W) derived from A is binary if 



W = -a ± + b ± (=|ox|+&l). 



(2.6) 



and the two types of intervals of the binary QL have the lengths a and b. Eq. (2^) prescribes how the inflated versions 
of the two intervals are divided into the original intervals. 
Let us introduce two positive parameters by 



£ := b/a, ( := -a±/6j 



(2.7) 



3 



Then, Eqs. ( |2.3| ) and (2.4) are rewritten as 

r(l £) = (1 £)M, 



(2.8) 



f(-C 1) = (-C 1)M. (2.9) 

Note that £ > 1 because a < b. Since M is an integer matrix, the row vector (1 £) must be also a left eigenvector 
correspon ding to the eigenvalue f of M. Hence we have £ = £ with £ := — l/£ being the dual to £ as defined in 



Appendix A 1 . It can be shown that the ratio of the frequency of occurrence of the interval a to that of (3 is equal to 
£. Although £ and £ belong both to the quadratic field Q[r], their dependences on r are different in general from the 
case of the Fibonacci lattice for which £ = £ = r (= tq). 

The projections of A onto E\\ and E± yield modules, A|| and Aj_, which are denseJD sets. They are isomorphic to 
A, and there are a natural bijection between any pair of the trinity, {A, A||, Axj.Ej A QL is nothing but a discrete 

subset of A||. The projection modules are equal to the Z-modules Z[£] and Z[£] (= Z[£]), respectively, where the 
physical space and the internal space arc scaled so that a = b± = 1, which we will sometimes assume. Remember 
that p + g£ € Z[£] and — p( + q £ Z[£] are combined by the bijection mentioned above. 

A well-known formula for a binary QL yields the position of the n-th site (lattice point) of Q(W) with W = 1 + £: 



+ (£-i) 



i + C 



(2-10) 



where <f> is the phase parameter depending on the position of the window in the internal space. 

It is important that the hyperscaling operation T transforms a canonical set, {a, b}, into another canonical set, 
{a', b'} := {Ta, Th}. The new canonical set yields another binary QL, Q(W'), where the new window W is related 
to the old one by W' — t~ 1 W . The sizes of the two types of intervals a' or (3 1 of Q(W) are given by a' := ra and 
b' := rb. The hyperscaling symmetry of A results in. that tQ(W) ~ Q(W), which means that Q(W) is a self-similar 
QL, whose ratio of the self-similarity is equal to t.£3 However, it is not obvious whether the two types of intervals a' 
and (3' are uniformly decorated by the two types of intervals a and (3 because the inflation matrix does not directly 
prescribe the order of the latter intervals in a' or (3' . The uniformity is evident for the case of the Fibonacci lattice 
as illustrated in Fig. ^. A proof of the uniformity for a generic case will be given in the next section. Note that the 
ratio of the self-similarit y (i .ft^ the scaling factor r) of a binary QL is restricted to a quadratic irrational of the form 
( |A3| ) given in Appendix ^jjjl 

Since the hyperscaling symmetry of the mother lattice is essential for the self-similarity of the relevant binary QLs, 
we shall confine our considerations to mother lattice with hyperscaling symmetries. Then, every binary QL obtained 
by the projection method from such a 2D lattice is self-similar, so that binary QLs appearing in the present paper 
are all self-similar. Therefore, we may call a self-similar binary QL simply as a binary QL. 

It is evident that T n with v n S Z transforms a canonical set of basis vectors into another. We say that two canonical 
sets are scaling-equivalent if one of them is transformed by T n with 3 n G Z into the other. Then, the relevant two 
binary QLs are mutually similar, and we need not distinguish them in a classification of QLs. The canonical windows 
of two QLs which are scaling-equivalent are related to each other by W' — r~ n W with 3 neZ. 

Let W be the set of all the canonical windows. Then, it is a discrete subset of the half line R + . It has the scaling 
symmetry, t^W = W, and is divided into scaling-equivalent classes. Two binary QLs with windows belonging to 
different scaling-equivalent classes are not mutually similar. 

The parameter £ characterizing a binary QL satisfies £ > 1 and £ > 0, so that is a quasi-reduced quadratic irrational 



(QRQI) as defined in Refs. |10|Ji]| and discussed in Appendix A 6. Conversely, the companion matrix of every QRQI is 
a quasi-normal UFM, which determines a 2D lattice with a hyperscaling symmetry and the relevant binary QL. Thus, 
there exists a bijection (one-to-one correspondence) between any pair of the three sets: i) the set of all the QRQIs, 
ii) the set of all the quasi-normal UFMs, and iii) the set of all the binary QLs. In what follows, a binary QL specified 
by a QRQI £ will be denoted as Q{£}. Since an inflation matrix is a quasi-normal UFM, it determines uniquely a 
binary QL (exactly, an LI class of binary QLs). Note that every quasi-normal UFM can be an inflation matrix of a 
binary QL. 



III. CLASSIFICATION OF BINARY SELF-SIMILAR QLS INTO MLD CLASSES 

If, between two QLs Q and Q', there exists a uniform local-rule for the transformation of Q into Q', Q' is said 
to be locally derivable fcom Q; moreover, if Q is also locally derivable from Q', we say that Q and Q' are mutually 
locally-derivable (MLD)H 
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According to a general theory of MLD relationship developed in Ref. |[ a necessary condition for two QLs to 
be MLD from each other is that they are derived by the projection method from a common mother lattice, where 
two mother lattices which diff er ir t-t-he scales of the physical space and the internal space are not distinguished as 
mentioned at the end of Sec. [I A O Therefore, two binary QLs which are MLD but are not mutually similar are 
derived by two canonical sets, {a, b} and {a', b'}, which are inequivalent with respect to the hyperscaling. Let X 
be a unimodular matrix which combines two canonical sets: 

(a' b') = (a b)X. (3.1) 



Then, the argument in Appendix A 4 results in that the ratio of the two intervals, £' := b' /a' , is related to the ratio, 
£ := b/a, by the modular transformation: £' = A _1 (£); the two ratios are modular equivalent. This is a necessary 
condition for two binary QLs, Q{£} and Q{£,'}, to be MLD from each other. The next problem is to show that this 
is also a sufficient condition. 



Using Eqs. ( |2.6| ) and Eq. (3.1), we can calculate the size of the canonical window W as W = \(t — u)( — v + w\ 
so that the difference W — W with W = 1 + C belongs to the module Z[£], where the relevant UFM X is assumed 



to be given by Eqs. ( A18). Hence_we can conclude by a general theory of the MLD-relationship among different QLs 
that and Q{£'} are MLD.L3 This completes a proof of the sufficiency of the condition above. We will show 

in the next section directly the MLD-relationship between and Q{C} by deriving an explicit substitution rules 

(SR) combining the two QLs. At all events, there exists a bijection between the set of all the MLD classes of binary 
QLs and the set of all the modular equivalent classes (MECs) of QRQIs. Therefore, we can specify an MLD class 
by the corresponding MEC, ((ko, hi, k n -i)), of QRQIs, or, equivalently, by the corresponding MEC of the 

quasi-normal UFMs, where n is the rank of the class. Thus, there exist N (= ko + k\ + ■ ■ ■ + fc n _i) binary QLs in the 
relevant MLD class, and the N binary QLs are: Q{£,^} with j = 0, 1, 2, • • • , n — 1 and k = 0, 1, 2, ■ • • , kj — 1, 
where the £^ are QRQIs defined in Appendix |A 6|. We shall call N the order of the MLD class. Other important 



parameters of the MLD class is m and e = (—1)™, where m is the common trace of the relevant UFMs. 

It is evident that each MLD class has its own mother lattice. That is, there exists a bijection between the set of all 
the MLD classes of binary QLs and the set of all the mother lattices with hyperscaling symmetries. 

The MLD-relationship between two QLs, Q{£,} and Q{£'}, in an MLD class is particularly simple for the case 
where the modular transformation, £' = A _1 (^), coincides with either of the two elementary modular transformations 



presented in Appendix |A 5|. The matrix U given by Eq. (All) is identical to the matrix M given by Eq. (2.2), and 



the two canonical sets in this case are related to each other by a' = b and b' = a + b. Then, the SR relating the two 



sets of the intervals is the same as the SR (2.1) of the Fibonacci lattice: 

a' = /3, = a/3. (3.2) 
This SR is closely related to the matrix U , so that we may call U the substitution matrix of this SR. On the other 



hand, the two canonical sets are related, in the case of the matrix S given by Eq. (All), to each other by a' = a and 
b' = a + b. Then, we can show by a similar geometrical argument to the one by which the SR above is derived that 
the SR for this case is given by 

a' = a, = f3a, (3.3) 
which is illustrated in Fig. 0. The substitution matrix for this SR is S. 



A simple geometrical rule for the choice between the two SRs (3.2) and (3.3) is given as follows: the rule (3.2) or 



(3.3) should be chosen according as {b, c} or {a, c} with c :— a + b is a canonical set; c is the diagonal vector 
of the relevant parallelogram. This provides us with an elementary procedure which allows successive generation of 
canonical sets of basis vectors. This successive procedure can be regarded, alternatively, as a recursive procedure of 
generating a series of "canonical" basis vectors, a^, i = 0, 1, 2, • • •, with the initial conditions, ao = a and ai = b. We 
shall call it an additive algorithm. The inverse procedure of it is a subtractive algorithm: one and only one of {d, a} 
and {d, b} with d := b — a is a canonical set, where d is the second diagonal vector of the relevant parallelogram. 
The additive algorithm (or the subtractive algorithm) may be called the diagonal (or anti-diagonal) algorithm. At 
any rate, we can generate a both-infinite set of canonical basis vectors T := {a^ | i £ Z}, and the i-th canonical set 
of basis vectors is given by {a^, a^}, where k is the maximum number under the condition that it is smaller than i 
and the sign of the second component of a& is opposite to that of a^. Hence there is a bijection between T and the 
set of all the canonical sets of basis vectors, and we can identify them. The hyper-scaling symmetry of A results in 
a.j + jv = T&i, so that T is divided into scaling-equivalent classes (series), which form an A-cycle, i.e., a cyclically 
ordered set of N elements. 

The properties of the trinity, {A, A||, A^}, is naturally succeeded by the trinity, {r, rii, T^}, where r n :— {ai \ i € 
Z} and r_L := {a^~ | i E Z} are the projections of T onto Eh and E±, respectively. In particular, L|| and T± have 
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the scaling symmetry, a,i + N — i~a,i and aj- +N = faj~, and are divided into N scaling-equivalent classes, which form 
iV-cycles. The former of the two is only composed of positive numbers but the latter is not. If the signs of the members 
are ignored, T± coincides with W because each member of T is the "anti-diagonal" of the parallelogram associated 
with a canonical set of basis vectors. If {a^, a^} is a canonical set above, we obtain £j = ai/ak and Ci = —a>k~/ a i~i 
which are periodic: £,:+at = £i and Ci+N — Ci- A 2D matrix is defined with the canonical set by Ai := (a& a^). 
It satisfies Ai + i — AiS if af- and aj- +1 have a common sign but Ai + i = AiU otherwise. It satisfies, furthermore, 
TAi = A i+N = AiM, where M is the companion UFM of £. 
We may write Ai = (a h)Xi, where 

x < - ( l :i ) < s - 4 » 

is a unimodular matrix. It is a quasi- normal UFM if i > 2 but the reciprocal of a quasi- normal UFM if i < 0, 
while Xi = I. It satisfies -Xj+jv = Xj,M. Since Ai — (a^ a^), we obtain a^ = «,a + Wib, ai = Ui + Wi£ G Z[£], 
and af- — — UiC + G Z[£]. The series of rationals, Wi/ui (or —U-i/iV-i), i = 1, 2, ■ ■ •, is just the series of best 
approximants to £ (or £) . The procedure to obtain approximants to an irrational on the basis of the geometry of a 
lattice has been established in the classical number theory.Efl 

It is sometimes convenient to change a single suffix numbering, Xi, into a double suffix numbering, x)i , with 
i = j + nN, where j = 1, 2, ■ ■ ■ , N and n G Z. Then a relation combining £i and iCj+jv turns to that combining 
Xn and x%\ x . Using the equality Xn^ = X^M n together with the Cayley-Hamilton theorem, M 2 — mM + el = 0, 
we can conclude that the quantities an , A„ , Xn \ (M\ and an satisfy two-term recursion relations which are 
isomorphic to x^l_i — mxn ~ e%n}-i', the solution of the recursion relation is written as_x n — F n x\ — eF n -iXo with 
F n being the generalized Fibonacci numbers associated with the quadratic irrationals rO Remember, however, that 
the double suffix numbering must not apply to £ and (. 

The N QRQIs in the MEC introduced above form an iV-cycle of modular transformations, and the relation between 
successive two members of the cycle is given by 



= ( 3.5) 



$ ) =U(tt 1) ) (3-6) 



if k > but 



with k = kj-i — 1. It follows that successive two members of the cycle of N binary QLs are related by the SR ( |3.2| ) 
or (3.3). Moreover, if two members of the iV-cycle are chosen arbitrarily, each of the two is derived from the other by 
a successive operation of the two types of the elementary SRs. That is, the two QLs are combined by composite SRs. 

It is appropriate at this point to define exactly the SR which combines two QLs. We shall distinguish for a while 
the two types of intervals a' or 0' from two types of sequences a a or ap composed of the two types of intervals a or 
[3. An SR is represented as a 1 = a a , /3' = crp. Let us denote a QL composed of the two types of intervals a or (3 as 
Q[a, /?]. If Q[a, (3] is represented as a sequence of two types of sequences a a or dp, there exists a new QL Q'[a', f3'], 
so that Q[a, (3] = Q'[a a , dp]. If the SR is used to derive Q'[a', (3'} from Q[a, 0\, we may call it a passive SR. 
Contrary to it, if it is used to derive the latter from the former, we may call it an active SR. The passive SR and the 
active one are the inverse procedure of each other. 

The substitution matrix of any SR is naturally a quasi- normal UFM. If the QL Q{£} is transformed by a passive SR 
into Q{£'}, the relevant two quadratic irrationals is related by the modular transformation as £ = X(£') with X being 
the relevant substitution matrix. If an SR with the substitution matrix X is operated first, and the second one with 
X' is operated subsequently, the substitution matrix of the resulting composite SR is given by XX' in the passive 
case but X'X in the active case. Remember that, for the passive case, the composition of two SRs is anti-isomorphic 
to the product of the two relevant substitution matrices. In the present paper, we will use SRs in the passive meaning 
unless stated otherwise. 

Each of the N scaling-equivalent classes in W has one and only one member in the fundamental interval, F := ]1, t\. 
Let Wo := YVHF with F := [1, r]. Then, it includes N+l windows. The pattern of the distribution of the windows in 
Wq characterizes completely the relevant MLD class, so that Wo is like a finger print of the MLD class. For example, 
if £ = (1, 2, 1) = (1 + Vl0)/3, we obtain W = {1, £ 1 + f, 1 + 2£, 2 + 3£}, which is shown in Fig. |. The last 
member of this "finger print" is equal to t = 3 + VTO , which is scaling-equivalent to the first member. 

In general, F is divided into N subintervals by Wo- Every subinterval is limited by two canonical windows, and the 
two relevant binary QLs have one common interval. If a window is chosen to be an internal point of the subinterval, 
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the resulting QL is a ternary QL composed of the three types of intervals associated with the two binary QLs. liilThat 
is, the ternary QL is a hybrid of the two binary QLs. 

The canonical windows are vertical widths of the fundamental parallelograms, and we may call their horizontal 
widths canonical interval (lengths) by a reason to be understood shortly. Let us denote by V the set of all the 
canonical intervals. Then it is identical to r || . The sizes of the two intervals included in a binary QL derived from 
the relevant mother lattice are given by consecutive two members of V. A new "finger print" Vo can be defined with 
respect to the scaling-invariant set V in a similar way as the one for Wo- It is remarkable that the two types of "finger 
prints" coincide for the MLD classes like ((1, 2, 1)). The reason will be revealed in a later section. 



IV. COMPOSITE SUBSTITUTION RULES AND INFLATION RULES 



It is sometimes convenient to employ composite SRs associated with the two types of UFMs which are given as 
blocked forms, L = S and K = US k ~ 1 with k and I being natural numbers. The SR for the former case is I fold 



repetition of the rule (3J3), so that we obtain 



= (3a l 



(4.1) 



where a 1 stands for a concatenation of I as. This SR is a direct conseque nce o f the equation s, a' = a and b' = /a + b. 
The explicit form of the substitution matrix of this SR is given by Eq. (A13) in Appendix A 2. On the other hand, 
the SR for the substitution matrix K = US^ 1 is a composition of the rule (H-l[ ) with I = k — 1 and the rule (|3.2]), so 
that we obtain 



(3, (3' = a(3 k . 



(4.2) 



This SR is a direct consequence of the equations, a' = b and b' = a + fcb. The two block SRs just obtained have 
s impl e geometrical meanings, which can also be understood easily. Since a quasi-normal UFM is repr esen ted a s Eq . 
( |A1C| ), every composite SR can be represented as a composition of a number of block SRs of the types (LI) and flt.St) . 

A composition of two SRs as well as a product of the corresponding two substitution matrices is incommutable. 
For example, the SR with the substitution matrix SU is given by 



t' = (3a, (3' — a(3a, 



(4.3) 



which is markedly different from the SR with the substitution matrix US; this SR is given by Eq. (4.2) with k = 2. 

Let Q and Q' be two QLs which are combined by the SR (4.3). Then, they are combined also by another SR, 
a' = a(3, (3' = a 2 /3, which is obtained from the former by moving the interval a at the right ends of the two relevant 
sequences of the intervals to their left ends. Therefore, these two SRs are equivalent. By similar procedures, we can 
obtain another two equivalent SRs: a' = (3a, (3 1 = (3a 2 and a 1 — a/3, (3 1 — a(3a. Conversely, by applying repeatedly 
the inverse ramcedure to the last SR, we can retrieve other three SRs. These proper and inverse procedures are called 
cyclic shifts.li-3 They yield a linearly ordered finite set of equivalent SRs, and the set is retrieved from any member of 
the set. The SRs in the set have a common substitution matrix. We will not distinguish equivalent but apparently 
different SRs. Note that the full set of equivalent SRs has a sort of mirror symmetry; an SR in the ordered set and 
the one at the mirror site are mutually mirror symmetric because each sequence of the former is the mirror image of 
the relevant sequence of the latter. 

An MLD class of binary QLs has the structure of an TV-cycle as mentioned in the preceding section. Therefore, if 
the two types of elementary SRs are applied N times in an appropriate order to an arbitrarily chosen member of the 
iV-cycle, the chosen member is retrieved. In other words, every binary QL is transformed into itself by a composite 
SR. Therefore, it is self-similar, and the composite SR is nothing but the IR. Since the inflation matrix of a QL is 



written as Eq. (A16), the relevant IR is obtained by a composition of the two types of the block SRs (4.1) and ( 



V. EXAMPLES 



A. MLD classes of rank one 



MLD classes of rank one form a series specified by natural numbers, and most of important ID QLs including the 
Fibonacci lattice belong to some MLD classes of this series. The MLD class is specified by ((m)) in our notation, for 
which N — m and e = — 1. It has only one RQI as given by 
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- (m + \J m 2 + 4^ , 



(5.1) 



i.e., the precious mean, which presents the scale of the self-similarity of the binary QLs in the class. There exist to 
QRQIs, £j: 0) = r-k with k = 0, 1, ■ • • , m - 1. 

The inflation matrix for is written with Eq. ( A16| ) as 



m— fc— 1 



fc fc(m — fc) + 1 
1 to — fc 



(5.2) 



The IR for Q{^°' ) } is given as a composition of two block SRs corresponding to US 1 



i-k-i 



and (42), we obtain 



and S k . Using Eqs. (4.1) 



(3a k , (3' = a((3a k ) 



k\m — k 



(5.3) 



which is equivalent to the IR given in Ref. |TJ or given by equation (20) in Ref. |[ 

There exists only one MLD class for AT = 1 or N = 2 but there exist two or more for N > 3. The MLD class ((1)) 
is only composed of the Fibonacci lattice, while ((2)) is composed of two MLD classes whose ratio of self-similarity is 
the silver mean, r = 1 + \/2 . The IRs of the two binary QLs in ((2)) are given by the above equations with fc = 
and fc = 1, and the relevant quadratic irrationals aref^ ' = 1 + V2 and ^ 1 °' > = \/2 , respectively. The two QLs are 



transformed to each other by the elementary SRs, (3.2) and (3.3) 



The ratio of self-similarity, r, of a QL must be generally a positive root of the quadratic equation z — raz + e = 
with m being a natural number and e = ±1. There exists at least one MLD class with given to and e. A representative 
of such MLD classes is ((to)) if e = — 1. There exists only one MLD class with given to and e if to < 5 and e = — 1. 
The two MLD classes, ((6)) and ((1, 2, 1)), are commonly characterized by to — 6 and e = — 1 (namely, r = 3 + vTO ). 

A remarkable feature of binary QLs of rank one is that the relevant parameters satisfy the condition, Z[£] = Z[r], 
which results in that the relevant QLs have simple IRs. 



B. MLD classes of rank two 

MLD classes of rank two form a 2D "series", ((fco, ki)), each of which is specified by a pair of natural numbers, 
and we obtain e = +1. We can restrict our considerations to the case fco > fci. An important subseries is formed of 
((to — 2, 1)) with m > 4, for which N = m — 1. This subseries has been investigated by several authors (see, for 
example, Ref. |l4| ). 

The companion matrix Mo — K\K^ of the RQI da := (to — 2, 1) is given by 

/ 1 to — 2 \ , r 

M °=(lm-lJ< ( 5 - 4 ) 

whose Frobenius eigenvalue is 

t = i (m + \/m 2 - 4 ^ , (5.5) 

i.e., the "anti-precious mean".i It presents the scale of the self-similarity of the binary QLs in the relevant class, and (9o 
is written as 9q = t — 1. The relevant MEC of QRQIs has m — 1 members, ^°' ) := t — fc with k = 1, 2, • • • , m— 2 and 
^q 1 "* = (1, to— 2) = (t— l)/(m— 2), where the parameter fc is shifted for a convenience from that following the convention 
in Sec. [T| The inflation matrices for the QLs, QU£ 0) } and Q{^ 1] }, are S ,fc - 1 C/ 2 S" n ~ fe - 2 and US m - 3 U {= 'Mo), 
respectively. The IR for Q{£ fc °^} is shown to be equivalent to 

a' = a k l3, ft = a k - 1 /3(a k P) m - k - 1 , (5.6) 

which is also equivalent to the IR given in Ref. |l4| or given by equation (20) in Ref. ||. On the other hand, the IR for 
Q{£o^} ^ s shown to be equivalent to 

a' = p m - 2 a , = f3 m ~ 1 a. (5.7) 
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If the notations for the two types o f int ervals are excha nged ia. this IR, the result coincides with the IR ( |5.6| ) but with 
k = m — 1. Therefore, the formula (5.7) is absorbed in (|5.6| ).E-J Note, however, that the second interval is shorter than 
the first one after the exchange of the notations. The explicit form of the inflation matrices are: 



M k = S^U'S™-*- 2 = ( \ k{m ~y~ 1 ) • (5.8) 

The simplest MLD class of the subseries ((m — 2, 1)) is given by setting m = 4. It includes three binary QLs whose 
ratio of self-similarity is t = 2 + ^/2>. 

The inflation matrices of binary QLs of rank two have the property, Z[£] = Z[r] or Z[£ _1 ] = Z[t], only if they 
belong to the subseries above but this property is not possessed by other binary QLs of rank two nor those of higher 
ranks. 

A representative of MLD classes with given m but with e = +1 is ((m — 2, 1)). There exist no other MLD classes 
if m < 7. The two MLD classes, ((7, 1)) and ((3, 2)), are commonly characterized by m = 8 and e = +1 (namely, 
t = 4 + vl5). In fact, the second MLD class is the simplest among rank two MLD classes which do not belong to 
the subseries above. It has five binary QLs corresponding to five QRQIs: 

l + VTS 
2 



The MLD class ((3, 2)) is considered to be the simplest member of the series of rank two MLD classes, ((k, 2)) 
with k > 3. The relevant parameters are: m = 2k + 2, e = +1, and N — k + 2. An important QL in the MLD class 

{(k, 2)) is <3{^ 1} }, whose inflation matrix is S(US k - 1 US)S- 1 = SUS^U. Explicitly, the inflation matrix coincides 
with a quasi-normal UFM defined by 

so that #o = (t — k — l)/fc. Since M is factorized into three blocks, (S)(US k ~ 1 ){U), the relevant IR is obtained as a 
composition of the three relevant SRs: 

a' = a{f3a) k , = a(/3a) fc+1 . (5.10) 

A series of binary quasiperiodic sequences with this type of IRs has been investigated in Ref. |lj. Note, however, that 
the sequences in Ref. O are not geometrical objects but those of letters. 



AO) _ 6 
so — — 

^ = 1 



£1 



(0) 



(i) 



i + vis 

2 ' 
15 



AO) 

?2 



C. MLD classes of rank three 



An important series of QLs formed of rank three MLD classes is given by ((1, k, 
parameters are: m — 2k + 2, e — — 1, and N = k + 2. The companion matrix Mq = 
9q := (1, k, 1) coincides with the quasi- normal UFM defined by 



Ml 



k k+1 
k + 1 k + 2 



1)) with k > 2. The relevant 
UKU = U 2 S k - x U of the RQI 



(5.11) 



so that 6*n 



(r — k)/(k + 1). The IR for Q{6 }, for example, is given by 

a' = {3(af3) k , (3' = /3(a(3) k+1 . 



(5.12) 



A series of binary quasiperiodic sequences with this type of IRs has been also investigated in Ref. [12]. The simplest 
MLD class in the series is ((1, 2, 1)), which has appeared at the end of the subsection VA. The inflation matrix 
given by Eq. (A2S) is a special case of the above with k — 2. 
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VI. SYMMETRIC SUBSTITUTION RULES AND SYMMETRIC INFLATION RULES 



The composite SR ( |4.2| ) with k = 2 is equivalent to 

a' = (3, /3' = (3a(3, 



(6.1) 



whic h is explicitly symmetric, and we may say that the SR (^^) with k = 2 is symmetric. M ore generally, the SR 
(4.2) is symmetric if k is even but not if k is odd. A similar argument applies to the SR (4.1). Explicitly, we may 
write 



a. 







a 1 ' (3a 1 ' 



if / = 21' with V being a natural number. In particular, if I = 2, it reduces to 

a' = a, = a(3a. 



(6.2) 



(6.3) 



On the other hand , the two e lementary SRs (3J2) and (3J3) are asymmetric. Also the SR ( |4.3| ) is asymmetric, while 
the t wo I Rs, (5.1C) and (5.12), are symmetric irrespective of the value of k. We should remark that the two SRs, (6.1) 
and (6.3), have simple geometrical interpretations in the projection method. 

We shall consider here the condition for an SR to be symmetric. A sequence composed of two types of intervals, a 
and (3, cannot be symmetric if both the number of a in it and that of (3 are odd. Thus, the parities of these numbers 
are important. It can be readily shown that a necessary condition for an SR to be symmetric is that its substitution 
matrix is a quasi-jpermutation-matrix (QPM) as defined in Appendix A 3 . This condition is shown to b e a sufficient 
condition as wello This is consistent with above obse rvations. For example, the inflation matrix (5.2) is a QPM if 
and only if both m and k are even, so that the IR ( |5.3| ) is symmetric only when this condition is satisfied. Explicitly, 
the symmetric IR is equivalent to a 1 = a k (3a k and (3' = (a') m ~ k a(a') m '~ k ' with m' := m/2 and k' := k/2. A 
similar argument applies to the IR (|5 



It is shown in Appendix A 3 that M with M being an inflation matrix is always a QPM if m is even. Since M is 



associated with the double IR, the double IR is equivalent to a symmetric IR. Similarly, the triple IR is equivalent to 
a symmetric IR if m is odd. Note that the scal ing ratio of the double or triple IR is equal to r 2 or r 3 , respectively. 
For example, if the symmetric IR given by Eq. (|4.3|) is doubled, we obtain the symmetric IR: 



= a{(3a) 2 , (3' = a{(3af 



(6.4) 



whose ratio of inflation-symmetry is (1 + \[2 ) 2 . On the other hand, if the IR of the Fibonacci lattice is tripled, we 
obtain the symmetric IR: a' = (3a(3 and (3' = (3(a(3) 2 , whose ratio of inflation-symmetry is Tq. 

The central interval in the sequence a' of the symmetric IR (5.12) is a and that in (3' is (3 if k is odd but these 
correspondences are reversed if k is even. The relevant substitution matrix is congruent in modulo 2 with / for the 
former case but with J := (J J) for the latter. These observations are general properties of symmetric SRs and 
symmetric IRs. 

If one of two QLs is derived from the other by a symmetric SR, and, moreover, the converse is also true, we may 
say that they are symmetrically MLD (s-MLD). Since the s-MLD is stronger than the conventional MLD, an MLD 
class can be divided into two or more s-MLD classes. For example, the MLD class ((m)) with an even m necessarily 
bifurcates because it includes two types of QLs with different scalings with respect to symmetric. IRs. Note that a 
necessary condition for the two QLs to be s-MLD is that the relevant two windows is concentric .EUl 

A symmetric SR is called elementary if it cannot be represented as a composition of two symmetric SRs. Such 



an SR is associated with an elementary QPM defined in Appendix A3. In particular, the SRs (6.1) and (6.3) 



are elementary symmetric SRs. Each s-MLD class is a cycle, and its two successive me mbe rs are com bined by an 
elementary symmetric SR. There are two important series of elementary QPMs, namely, (5J3) and ( 5.11 ) with k > 0. 
This is confirmed by the expressions M+ = SUS k ^U and M T = U 2 S k - 1 U for k > 1 but A/ + = S' 2 and M ~" = US. 
The corresponding symmetric SRs are identical to the IRs, ( 5.10 ) and ( 5.12 ), respectively. That is, QLs belonging 
to the two series of QLs with inflation matrices, with k > 1 and Mr with k > 0, have the property that their 
inflation matrices are elementary. The two series are the main subjects in Ref. [l2[ 

From here to the end of this section, we use SRs in the active meaning for a convenience; the directions of the 
arrows to appear should be reversed in the case of the passive meaning. A simplest example is again the case of a 
rank one MLD class, ((ni)). For brevity we set Qk '■= Q{£,^}, whose inflation matrix is denoted as M^. If m is even, 
the inflation matrix for Q a is factorized into elementary QPMs as M (= US" 1 ^ 1 ) = (Mr^M^)™/ 2-1 , which yields 
the following cycle of symmetric SRs: 
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Qo -► Q 2 ~> ► Qm-2 Qo, (6.5) 

where the single arrow stands for the SR ( |6.3[) but the double arrow the SR (6.1). The above cycle only includes the 



QLs Qk with even k. For Q& with an odd k, the double IR must be employed. The relevant inflation matrix for Qi 
is written as (Mi) 2 = SUS'^US" 1 - 2 = (M+)(M^) m / 2 " 1 , so that we obtain the following cycle of odd k QLs: 

Qi -> Qi -* ► Qm-i => Qi, (6.6) 



where the double arrow stands for the symmetric SR which is given by ( 5.10 ) with k = to. Thus, the MLD class 
((to)) with an even to is divided evenly into the even s-MLD class and the odd s-MLD class. Each of the t wo s-MLD 
classes for the c ase of to = 2 is composed of a single QL, Qo or Qy. The symmetric IR for Qo is given by ( |6.l|) , while 



that for Qi by (6_J) 



We shall proceed to the case of the MLD class ((to)) with an odd to. Then triple IRs are in order. The relevant 
inflation matrix for the case of Qo is factorized as 

(M ) 3 = (t/S™- 1 ) 3 

= (M -)(M+)(™- 3 )/ 2 (M+)(M+)(— (6.7) 



which reduces to Eq. (A.17) if to = 5. Since the right hand side has to factors, the to QLs, Qk with k = 0, 1, • ■ • , to— 1, 



are cyclically related by elementary symmetric SRs: 

Qo Q2 —> ■ ■ ■ —> Qm-i Qi — * Q3 — * • • • 

-» Qm-2 => Qo- (6.8) 

That is, this MLD class is simultaneously an s-MLD class. 

The above results for the case of a rank one MLD class is readily extended to the case of the subseries ((to — 2, 1)) 
of rank two MLD classes. We only present here the results of factorizations of the relevant inflation matrices. 
The case with an even m has two s-MLD classes, and the relevant factorization for the even s-MLD class is M2 = 
(M+)(M+) m / 2 - 2 but the one for the odd s-MLD class is (Mi) 2 = (M-_ 2 )(M ')(M+) m ^ 2 - 2 The relevant factorization 
for an odd to is (Mi) 3 = (M m _ 2 )(M Q )(M+^ m ~ 5 ^ 2 (M+)(M+Y m -^/ 2 . Note that Q k with k = is forbidden in 
the present case, so that the even s-MLD class with even to is composed of to/2 — 1 QLs. The situations are more 
complicated for other MLD classes of rank two or of higher ranks. 

A quick derivation of a symmetric IR is possible if we use blocked SRs. The SR ( |6.2| ) is associated with the block 
(S 2 ) 1 ' . Another two SRs associated with the two important blocks (M / t)(M f ) p and (M A T)(M + ) P are written with 
$ := aPpa p as 

a' = p(afj) k , (3' = (3(a/3) k+1 , (6.9) 



*' = a(/3a) k , (3' = a0a) 



k+l 



(6.10) 



Although the SR (4.3) is asymmetric, it is equivalent to the symmetric but fractional SR: 



(6.11) 



It can be shown generally that every asymmetric SR is equivalent to a symmetric but fractional SR.Q The same is 
true for asymmetric IRs. Thus, an MLD class is simultaneously an s-MLD class if fractional SRs are allowed. 



VII. MISCELLANEOUS SUBJECTS 



A. Periodic approximants 

Let a 1 and j3' be two types of sequences of intervals a and (3 of a QL, Q, and assume that the pair of the sequences 
forms an SR. Then, the periodic structure (a 1 ) 00 (or (/3')°°) obtained by an infinite concatenation of a' (or j3') is 
called a periodic approximant to QO It is evident that the period of a periodic approximant is equal to a canonical 
interval, and different periodic approximants have different periods. Therefore, we shall denote by Qn (n > 0) the 
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periodic approximant corresponding to the canonical interval a„ . The periodic approximant Q„ is changed by an 
active IR to Q„l_i- Hence all the periodic approximants are derived by the successive application of the IR to the 

prototypes, Q^f 1 with j = 0, 1, 2, • • • , N — 1.0 The first three prototypes are a, /?, and a/3, whose periods are 1, 
£, and 1 + £, respectively. The set of the periods of N prototypes is identical to the "finger print" Vq. There is a 
bijectioe, between the set of all the periodic approximants and the set of all the best approximants to the relevant 
QRQI.B 

So far we have confined our argument to periodic approximants of a single QL. If an SR is applied to a periodic 
approximant to the QL, we obtain a periodic approximant to another QL in the same MLD class to which the original 
QL belongs. 



B. Duality 

The properties of a QL are decisively ruled by the relevant quadratic irrational £. A quadratic irrational has the 
sole algebraic conjugate, and to take the algebraic conjugate of a quadratic irrational in the quadratic field Q[£] is an 
automorphism of Q[£]. This automorphism is a dual relationship, so that the QL may have some properties associated 
with the duality. 

If the mother lattice A is a square lattice, it has a canonical set of basis vectors {a, b} satisfying the conditions: 

a 2 = b 2 , a-b = 0. (7.1) 



If a and b are represented as Eq. (2.5), this condition leads aa± + bb± = 0, which is equivalent to £ = £ because of 



Eq. (2.7). Since £ = £, this means that £ is strongly self-dual. Conversely, if £ is strongly self-dual, the condition (7.1) 
is satisfied provided that the scale of the internal space is chosen appropriately. It follows that the mother lattice of 
an MLD class of binary QLs is taken to be a square lattice if and only if the relevant MEC of the normal UFMs is 
strongly self-dual. 

If £ is an RQI, so is its dual £. We shall denote the mother lattices associated with the two RQIs as A and A, 
respectively; the two are dual to each other. Also, W and V (or Wo and Vo) are dual to each other. Since the 
projection modules of the two mother lattices are given by Ay = Z[£] and A = Z[£], the two mother lattices coincide 
if and only if the relevant MEC of the normal UFMs is self-dual; we may say then that A is self-dual. If A is self-dual, 
we obtain W = V and Wo = Vo- We will discuss below the role of the dual mother lattice A. 

Let V be a finite interval in the physical space, and let Q(V, W) := Q(W) n V. Then, Q(V, W) is a projection of 
A n R onto the physical space, where R :— V x W, the direct product, is the rectangle obtained as a cut of the strip 
used to derive Q(W) in the projection method; we may write Q(W) — Q(R, W) with R being the entire physical 
space. Let Q(V, W) be the projection of An R onto the internal space. Then, we have Q(V, W) = Q(V) D W with 
Q(V) := Q(V, R), where R is the entire internal space. It is evident that Q(V) is a ID QL in the internal space, and 
the mother lattice of Q(V) is A. We may say that Q(W) and Q(V) (or, also, Q(V, W) and Q(V, W)) are dual to 
each other. 

The QL Q(V) is binary, if and only if V is a canonical interval. The two QLs in the dual pair can be mutually 
similar if and only if A is self-dual. At any rate, the dual relation between RQIs is associated with a symmetry 
between the physical space and the internal space. Note that Q(V, W) and Q(V, W) have a common number of 
lattice points. Moreover, they can be building blocks of the relevant periodic approximants. The present discussion 
on the dual relation applies to QLs in higher dimensions as well. The mother lattices of octagonal, decagonal, and 
dodecagonal QLs in 2D are shown to be all self-dual and so are for the cases of the three types (i.e., P, I, and F) of 
icosahedral QLs in 3D. 

The segments Q(V, W) with different canonical pairs {V, W} are not necessarily independent because the 
hyper-scaling symmetry of A yields Q(V, W) = T n Q(r~ n V, r n W) for v n E Z. Similarly we have Q(V, W) — 
T n Q(r n V, r~ n W). Note that the area of the rectangle, R := V x W, is invariant against the hyper-scaling. 

The structure factor of a QL, Q, is derived from the Fourier transform Q* of the QL, where Q and Q* are treated 
as Dirac measures (generalized functions) .113 Let A* be the reciprocal lattice to A. Then, the support of Q* is the 
Fourier module, which is the projection of A* onto the reciprocal physical space E*y Let A :— (a b) be a 2D matrix 
formed of the two basis vectors of A and A* := (a* b*) the counterpart in the reciprocal space. Then, we may write 

A= (-cQ' A * =a *(-£i) W 

with a* := 2n/(l + ££) because ( t A*)A = 2ttI. Thus, A and A* are dual to each other except a scale factor. The 
effect of the hyper-scaling on to A* is represented as TA* = A* t M. The structure factor of Q is closely related to a 
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ID QL whose mother lattice is A*. As a consequence, it is approximately self-similar. If A is self-dual in particular, 
the structure factor has a similar strpcture to the QL itself. Note, however, that A and A* are not dual for the case 
of icosahedral QLs of types I and 



VIII. SUMMARY AND DISCUSSIONS 



We have investigated ID binary QLs with geometrical self-similarities, and found that there exists a bijection 
between any pair of the four sets: i) the set of all the MLD classes of binary QLs, ii) the set of all the mother lattices 
with hyperscaling symmetries, iii) the set of all the MECs of quasi-normal UFMs, and vi) the set of all the MECs of 
QRQIs. If a mother lattice with a hyperscaling symmetry is fixed, the following the six A^-cycles are isomorphic to 
one another: i) the relevant set of all the MLD classes of binary QLs, ii) the relevant MEC of quasi-normal UFMs, 
iii) the relevant MEC of QRQIs, iv) the set of all the scaling-equivalent class of the canonical basis vectors, v) a 
similar set of the canonical windows, and vi) a similar set of the canonical intervals. The algorithm which determines 
the structure of the cycle has been independently given for the three cycles, ii), iii), and iv), although the three 
algorithms are equivalent. The algorithm which determines the structure of the cycle iii) is identical to the algorithm 
for the continued-fraction expansion of a QRQI, £■ It is equivalent to the anti-diagonal algorithm which determines 
the structure of the cycle iv), while the diagonal algorithm is identical to the algorithm for the continued-fraction 
expansion of a QRQI, £. Two successive members of every MLD class of binary QLs are combined by one of the two 
elementary SRs, and arbitrary chosen two members are combined by composite SRs. The IR of a QL is given as a 
special composite SR. 

Every ID binary QL has a symmetric IR. A necessary and sufficient condition for an SR to be symmetric is that 
its substitution matrix is a quasi-permutation-matrix (QPM). A conventional MLD class can be divided into two or 
more s-MLD classes. There exists a bijection between the set of all the s-MLD classes of binary QLs and the set of 
all the MECs of QPMs. Each s-MLD class is a cycle, and its two successive members are combined by elementary 
symmetric SRs. 

There exist two quadratic irrationals, £ and r, characterizing an MLD class of binary QLs. Of the two, £ includes 
complete information on the MLD class but r does not in general, so that the former is more fundamental. The latter 
has complete information only when the case of MLD classes of rank one or of a special series of MLD classes of rank 
two. 

The binary QLs investigated in the present paper have IRs whose inflation matri ces are unimodular. Moreover, 
the IRs are compositions of the elementary SRs of the two types given by Eqs. (3J3) and (3.2). Then, general 
arguments given in Refs. || and ^ result in that the relevant trace maps have the invariant of the standard form: 
/ := x 2 + y 2 + z 2 — 2xyx — 1. Therefore, the energy spectra and the one-electron wave functions of the tight- 
binding models on these structures are fractals which belong to the same classes as those of the Fibonacci lattice.EJ. 
Furthermore, we can conclude that each MLD class of binary QLs has its own universality class with respect to the 
one-electron properties because two successive members of every MLD class are combined by one of the two elementary 
SRs.il. 

We have shown that every MLD class is coded by a cycle of natural numbers. However, the MLD class is an N- 
cycle, and is not isomorphic to the cycle coding the MLD class. Therefore, we shall pursue a more convenient coding 
system. Let us take the MLD class ((3, 2, 4)) as an example, and transform it into the binary code 100010100, which 
is a reversed concatenation of the three numbers 100, 10, and 1000 with three, two, and four digits, respectively. If 
the digit 1 in the binary code is replaced by the elementary UFM U and the digit by S, we obtain a quasi-normal 
(exactly normal) UFM which is the inflation matrix of a member of the MLD class. The inflation matrices of other 
members of the MLD class are obtained similarly from cyclic permutations of the binary code. Therefore, for every 
binary QL, we shall associate a cyclic binary decimal in the interval ]0, 1] so that the first period codes the relevant 
inflation matrix. Since we can naturally define a cyclic equivalence between two cyclic decimals, we can conclude that 
i) there exists a bijection between the set of all the binary QLs and the set of all the cyclic binary decimals in the 
interval ]0, 1], ii) there exists a bijection between the set of all the MLD classes and the set of all the cyclic equivalence 
classes of cyclic binary decimals, and iii) the two cycles related in this bijection is isomorphic. Thus, an MLD class 
of order N is coded by an irreducible cycle of N binary digits like ((100010100))b, where the suffix b stands for the 
symbol being a binary code. 
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APPENDIX A: MATHEMATICAL GLOSSARY 



M ost r esul ts to be presented in this appendix are found in the literature, e.g., Ref. Ea but several results in appendices 



A 2 , A 3 and A 6 seem new 



1. Unimodular matrix 

A matrix M whose matrix elements are all integers is called unimodular if det M = e = ±1. Let {a, b} a set of 
basis vectors of a 2D lattice and {a', b'} be another set. Then, the two sets are related to each other by a unimodular 
matrix 

«=(") (A1 » 

as 

(a' b') = (a b)M. (A2) 

Conversely, if M is unimodular, a' and b' form a set of basis vectors of the 2D lattice. 

There exists a linear transformation T so that a' := Ta and b' := Tb; T is an automorphism of the 2D lattice. 
Since the signs of basis vectors of a lattice are not very important, we can assume that m = TrM is a nonnegative 
integer. The eigenvalues of M are the roots of the quadratic equation z 2 — mz + e = 0. There can be three cases, 
elliptic, parabolic, and hyperbolic, depending on whether the discriminant D :— m 2 — 4e is negative, zero, or positive, 
respectively. We are interested in the hyperbolic case, to which a general argument to be given hereafter will be 
confined. Then, m must obey the condition m>life = — 1 but the condition m > 3 if e = +1. Of the two roots of 
the quadratic equation, the one which is larger than unity is given by the real quadratic irrational: 

t = - (m + \/m 2 - 4e) . (A3) 

The second root is given by f, the algebraic conjugate of r, and r and f satisfy 

tt = e, t + f = m, (A4) 

so that |f | < 1. An important Z-module is generated by r: 

Z[r] := {t + ur\t, ue Z}. (A5) 

Remember that the Cayley-Hamilton theorem yields M 2 — mM + el = with / being the unit matrix. 

The linear transformation T has two principal axes which we can assume as shown shortly to be mutually orthogonal; 
T enlarges along one of the axes by r but shrinks along the other by r" 1 = \f\. If the coordinate system matches the 
principal axes, T is represented by the diagonal matrix 

r-(;;). (A6) 

Let A :— (a b) be a 2D matrix formed of the two column vectors represented in this coordinate system. Then, Eq. 
( |A2| ) is equivalent to 

TA = AM, (A7) 

so that M is diagonalized by A: AMA^ 1 — T, and the two row vectors of A are the left eigenvectors corresponding 
to the two eigenvalues, r and f , respectively. Therefore, the basis vectors a and b are almost determined from M 
by the orthogonality of the principal axes; we have a freedom of choosing proportionality constants of the two left 
eigenvectors. _ . 

Let (a b) be the left eigenvector given by the first row of Alj Then, it satisfies the equation, 

r(a b) = (a b)M. (A8) 

Hence the ratio £ := b/a satisfies r(l £) = (1 £)M. It follows that 
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r = p + r£, r£ = q + s£. 



(A9) 



Hence £ is a real quadratic irrational belonging to the quadratic field Q[t]. It is important that the Z- module Z[£] 
coincides with Z[r] only when r = 1 because Z[t] = rZ[£]. Note that £ is a root of the quadratic equation: 

r£ 2 + (p-a)£-g = 0. (A10) 

In this paper, we always mean by a quadratic irrational a real quadratic irrational. 

The ratio £ is uniquely determined by M but the converse is not necessarily true because M n with n being any 
nonzero integer satisfies the condition. A unimodular matrix is called irreducible, if it cannot be written as a power 
of another unimodular matrix. Every quadratic irrational £ has its proper irreducible unimodular matrix; we shall 
say that £ and the relevant matrix are companions of each other. 

Another quadratic irrational is defined with £ by £ := — We shall call it the dual to £ because the tilde operation 
is recursive. The companion matrix of £ is equal to t M, the transposed matrix of M. 

The set, GL(2, Z), of all the 2D unimodular matrices form an infinite but discrete group. The set, GL'{2, Z), of 
all the hyperbolic and irreducible members of GL{2, Z) is, however, not a group because it does not include the unit 
matrix. We may say generally that two unimodular matrices M and M' are modular equivalent if there exists a 
unimodular matrix K so that M' = K~ X MK. The set GL'{2, Z) can be grouped into modular equivalent classes 
(MECs), which are disjoint. The two eigenvalues are common among different members of an MEC. Also the trace 
and the determinant are common, so that they are proper numbers of the class. Note, however, that the pair of the 
two numbers is not sufficient to specify an MEC. For example, the two matrices, (° g) and Q ^) are shown to be 
not modular equivalent. 

There exists a bijection (one-to-one correspondence) between any pair of the three sets: i) GL(2, Z), ii) the set of 
all the sets of basis vectors of a 2D lattice, and iii) the set of all the automorphisms of the 2D lattice. Therefore, 
we can define primed subsets like GL'(2,Z) for the latter two sets as well, and the division of GL'(2,Z) into MECs 
induces similar divisions of other two primed sets. 



2. Unimodular Frobenius matrices 



The case where M £ GL'(2, Z) is a Frobenius matrix, whose matrix elements are all nonnegative, is of particular 
importance. The positive eigenvalue, r, of M is called the Frobenius eigenvalue, and the relevant left or right 
eigenvector is called the Frobenius eigenvector, which is assumed to have two positive components. It follows that 
the companion quadratic irrational £ is positive. 

The set of all the Frobenius matrices in GL(2, Z) do not form a group but only do a semigroup because the inverse 
of a unimodular Frobenius matrix (UFM) has not necessarily the Frobenius property. Let M and M' be two UFMs 
which are modular equivalent to each other. Then, we say that they are strongly modular equivalent to each other 
if there exist two UFMs K and L so that M' = K~ X MK and M — L~ X M'L. This condition is equivalent to the 
one that M is a product of the two UFMs, M = KL. More generally, if M is a product of a number of UFMs, it is 
strongly modular equivalent to another UFM which is obtained from M by a cyclic permutation of its factors in the 
product. 



A UFM represented as Eq. (Al) is called quasi-normal if it satisfies the two conditions: i) it is irreducible and ii) 
p < q and r < s. It is called normal if it satisfies the condition, iii) p < r and q < s, in addition to the former two. A 
quasi-normal UFM is normal if and only if its transpose is also quasi-normal. In particular, a UFM and its transpose 
are simultaneously normal. Note that a normal UFM is always hyperbolic but a quasi-normal UFM can be parabolic. 

The set of all the quasi-normal (or normal) UFMs in GL(2, Z) form a semigroup if the unit matrix is included in 
it. If an irreducible UFM M /) is not quasi-normal, P~ 1 MP with P := (° ) is quasi-normal. That is, a UFM 
being not quasi-normal is strongly modular equivalent to a quasi-normal UFM. Therefore, properties of the former are 
investigated from those of the latter. An important property of a quasi-normal UFM is that its companion quadratic 



irrational is larger than unity, which follows from Eq. (A£). 
Simplest but nontrivial quasi-normal UFMs are: 



s=( i m, u= ? ; i. .aid 



1 J ' ^11 

In fact, the second of the two is normal. A scries of normal UFMs is defined by 



A H H) ,A,2, 
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with k being natural numbers. We shall call a member of the series an elementary normal UFM. The matrix ( A12j ) 
is factorized as K = US k ~ 1 . Similarly, a series of quasi-normal UFMs is defined by 

l =G0 (A13) 

with I being natural numbers. We may wr ite L = S l , which is a generic form of a parabolic quasi-normal UFM. 



If a normal UFM is represented as Eq. (Al), we obtain 



MR' 1 = ( P ) (A14) 

1 s — kr r ' 



We shall denote by fco the maximum number k under 
UFM. Then, MKq 1 is a normal UFM unless it is equal to the unit matrix. That is, a normal UFM is transformed 
into another normal UFM with a smaller "norm" provided that the norm is defined appropriately. We can repeat this 
procedure until we arrive at an elementary normal UFM. Thus, a normal UFM is uniquely factorized into elementary 
normal UFMs: 

M = K n _\K n _i ■ ■ ■ K Q . (A15) 

The rank of M is defined by the total number n of the factors in the right hand side of this expression. We should 
mention that the algorithm of this factorization is nothing but a parallel Euclidean algorithm. We can specify M by the 
ordered set (fco, k%, • • ■ , k n -\) of natural numbers, where fc, is a number specifying Kf, we have adopted the reversed 
order for a convenience of a later argument. The members of the ordered set are not entirely independent because of 
the irrcducibility condition i) above. For example, fco an d k\ must be different if n — 2. If M — (kg, k%, • • • , fc„_i), 
then M = (fc n _i, fe„_2, • • • , fco); M is symmetric if and only if the relevant ordered set is mirror-symmetric. 

As a result of the uniqueness of th e factorization above, we can conclude that a normal UFM which is modular 



equivalent to the normal UFM ( |A15| ) has a similar factorization to above but the factors are cyclically permuted. 
Hence, the rank is common between the two UFMs. An MEC of normal UFMs is composed of n members, where 
n is the common rank of them. The MEC is specified by the cyclically ordered set ((fco, fci, ■ ■ ■ , fc n _i)), where 
double angular brackets are used to distinguish the set from the linearly ordered set above. Note that the sign of the 
determinant of every member of the MEC is equal to (—1)™. We will call hereafter a cyclically ordered set simply a 
cycle. 

An MEC is called self- dual if it is specified by a cycle with a mirror symmetry. Every member of a self-dual MEC 
is modular equivalent to its transpose. An MEC is always self-dual if its rank is less than three. An MEC of rank 
three is self-dual if and only if it takes the form ((fc, fc', fc)), while an MEC of rank four is self-dual if and only if its 
form is one of the two alternatives, ((fc, fc', fc', fc)) and ((fc, fc', fc", fc')). 

An MEC is called strongly self-dual if it includes a symmetric member (i.e., a symmetric matrix). A self-dual MEC 
is strongly self-dual if its rank is odd but not necessary so if its rank is even. For example, an MEC of rank two is 
never strongly self-dual, while an MEC of rank four is strongly self-dual if and only if its form is the former of the 
two alternatives above. 



3. Quasi-normal UFM and quasi-permutation-matrix 



Let M' be a quasi-normal UFM. Then, we can show by a similar argument to that presented after Eq. (A14) that 
there exists the maximum number I under the condition that L~ X M' remains a UFM. In the maximal case, L~ Y M' 
is a normal UFM. Then, we can show that M := L~ 1 M'L is a normal UFM as well, so that a quasi-normal UFM 
is strongly modular equivalent to a normal UFM. It follows that a quasi-normal UFM is represented with a normal 
UFM M as S l MS~ l with I being a nonnegative integer. More precisely, if M is represented as Eq. (A15), we obtain 
fco quasi-normal UFMs: 

M k = S k US k — 1 - 1 US k — 2 - 1 ■ ■ ■ US kl - 1 US k °- k - 1 (A16) 

with k = 0, 1, 2, • • • , fco — 1 and Mq := M, where use has been made of the equations, K{ — U S ki ~ Y . Note 
that the factorization of Mj. like this is unique. More strongly, there exists an algorithm by which each factor in 
such a factorization is determined step by step: the right end factor of a quasi-normal UFM M is S if MS^ 1 is 
quasi-normal but U otherwise. The two elementary UFMs, U and S, are two of the four which are given in Ref. @ as 
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generators of UFMs. The remaining two are not necessary here because the UFMs under consideration are restricted 
to quasi-normal ones. 

If ((fco, ki, ■ ■ ■ , fc„_i)) is an MEC of normal UFMs, we can obtain fcj quasi-normal UFMs from the i-th member 
of the class; the total number of quasi-normal UFMs obtained is equal to N := fco + k\ + ■ ■ ■ + k n -\. They form an 
MEC of quasi-normal UFMs, which are strongly modular equivalent to one another. The two MECs, one normal and 
the other quasi-normal, can be specified by a common cycle of natural numbers although the l atter MEC is a cycle 
composed of N members. If a member of an MEC of quasi-normal UFMs is factorized as Eq. ( A16| ), n, the number 



of U, and N, the total number of the factors, are proper numbers of the MEC; n is nothing but the rank of the MEC. 
We shall call N the order of the MEC. 



The subject in the remaining part of this subsection concerns only in Sec. VI in the text. There are only two 
permutation matrices in 2D, i.e., the unit matrix / and J := (® Q ). We shall call a UFM a quasi-permutation-matrix 
(QPM) if it is congruent in modulo 2 with / or J. The product of two or more QPMs is also a QPM. If a UFM is not 
a QPM, it is congruent in modulo 2 with one of the four matrices, S, 'S 1 , U, and R:= ( 1 ) because detM = ±1. In 
what follows, the congruence relation x = y mod 2 is simply written as x = y. Then, the four matrices just introduced 
satisfy S 2 = QS) 2 = I,U 2 = R, R 2 = U, and U 3 = R 3 = I. 

Since m = or m = 1 according as m :— Tr M is even or odd, respectively, a UFM with an odd m cannot be a 
QPM. More precisely, M = U or M = R for this case, so that M 2 is not a QPM but M 3 is; M 3 = I. On the other 
hand, a UFM with an even m satisfies M = S or M = l S if M is not a QPM. Then, M 2 = I. 

A QPM is called elementary, if it cannot be represented as a product of two QPMs. Every QPM is uniquely factorized 
into elementary QPMs. The two QPMs, S 2 and US, are simplest elementary QPMs which are quasi-normal. We 
present other four types of elementary and quasi-normal QPMs: SUS k U = I, SUS l U = J, US l US l U = I, and 
US l US k U = J, where fc is an odd integer but I and V are even. There can be more complicated QPMs, e.g., 
U(US k ) k U with fc and fc' being any odd numbers, but we shall not pursuit them. A QPM is called irreducible, if it 
cannot be written as a power of another QPM. If a composite QPM is factorized into elementary QPMs, we obtain 
another QPM by a cyclic permutation of its factors, and the two QPMs are strongly modular equivalent to each other. 
Note that an irreducible QPM can be reducible as a usual UFM. For example, if M is not a QPM, M 2 or M 3 is an 
irreducible QPM. 

Finally, an example of a factorization of a QPM into elementary QPMs is presented: 

(US 4 ) 3 = {US)(S 2 )(SUS 4 U){S 2 ) 2 . (A17) 



4. Modular transformation 



Let us take a unimodular matrix 



X= I * " I g GL(2, Z). (A18) 



v w 



Then, a transformation from an irrational number £ to another irrational £' is uniquely defined by the condition 

(1 = A(l £')X or A(l O = (1 (A19) 

where A is a parameter determined by this condition. This transformation is called the modular transformation, which 
is written explicitly as 

= ^7-^- (A20) 

Our definition of the modular transformation is slightly different from the conventional one but is substantially 
equivalent to the latter. If £' = X(£) and £" = X'(g') with X' being the second unimodular matrix, we obtain 



£" = X"(£) with X" := X'X, which can be readily proved from Eq. (A19). Therefore, the set of all the modular 
transformations form a group which is homomorphic with GL(2, Z). More precisely, it is isomorphic with GL(2, Z)/Z2, 
where Z 2 := {/, —J} is a normal subgroup of GL(2, Z) because two unimodular matrices X and —X yield an identical 
modular transformation. 

A quadratic irrational is changed into another quadratic irrational by a modular transformation. Hence, the set of 
all the quadratic irrationals can be divided into MECs, which are disjoint. 
Let us introduce with the unimodular matrix X a new set of basis vectors by 
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(a' b') = (a b)X. (A21) 



The first row of this equation reads 

(a' b') = (a b)X. (A22) 

It follows that the ratio £' := b'/a' is related to £ by the modular transformation, £' = X^ 1 ^). As a consequence, we 
can conclude that £ is a fixed point of a modular transformation associated with its companion matrix M: £ = Af(£). 
Moreover, we obtain £' = M'(£') with M' := X~ l M X. Thus, there exists a bijection between the set of all the MECs 
of quadratic irrationals and the set of all the MECs of (hyperbolic) unimodular matrices. 

The following proposition provides us with a different view for modular equivalence between two irrationals, £ and 
£': a necessary and sufficient condition for £ and £' to be modular equivalent is that the two Z-modules, Z[ £] an d 
Z[£'], are scaling-equivalent: Z[f] = /uZ[£] for 3 /x e R; if f = X" 1 ^) with X being the unimodular matrix ( [A1S| ), 



the multiplier /i is determined as u-^= t + v£. The scaling equivalence between Z-modules is essential in a classification 
of quasilattices in any dimensions .E3 



5. Number theory of continued fractions 

Let £ be an irrational number being larger than unity. Then, we can derive a second irrational £' > 1 by the 
equation: 



(A23) 



where the symbol [*] stands for the maximum integer among those which do not exceed the number *. This equation 
is written with the use of the modular transformation as 



(A24) 



where K is the normal UFM ( A12). Since K = US k 1 with S and U being UFMs (All), this modular transformation 
is represented as a composition of the two types of elementary procedures: 



£' = S(0=£-1, ? = U@ 



1 



£-1 



(A25) 



The first of the two procedur es is used when £ > 2, while the second when 1 < £ < 2. 

If the recursive procedure ( A23) is repeated indefinitely, £ is expanded into an infinite continued- fraction. The i-th 
step of this recursive procedure is given by 



Ci — ki 



1 



with £o :— £ and ko ■= k, so that 



£,o = k + 



1 



1 



fci + k 2 



1 1 

-i + Ti 



(A26) 



(A27) 



This equa tion is equivalent to £j 



M^ ) with Mi = h A, , 



■ Ko being a normal UFM, where Kj is given by 
Eq. ( A12 ) with k = kj. If £, in the continued- fraction ( A27 ) is replaced byr-an integer k satisfying 1 < k < ki, we 
obtain a rational number, which is shown to be a "best" approximant to £j. Ii3 

Numbers having periodic continued- fraction expansions are important. Let M be the normal UFM derived from 
the first period of the expansion of such number, £. Then £ is the companion quadratic irrational of M. Hence, £ 
must be a quadratic irrational. More strongly, £ > 1 and £ > 1 because M is normal. However, the continued-fraction 
expansion of a generic quadratic irrational is not always periodic; exactly, it is periodic except for first several terms. 

A quadratic irrational is called a reduced quadratic irrational (RQI) if its continued-fraction expansion is periodic. 
A necessary and sufficient condition for a quadratic irrational 9 to be an RQI is that it satisfies the conditions: 9 > 1 
and 9 > 1, which is equivalent to — 1 < 9 < 0. Thus, there exists a bijection between the set of all the normal UFMs 
and the set of all the RQIs. Then, the symbols used to specify each normal UFM and each MEC of normal UFMs 
are used to specify the corresponding RQI and relevant MEC of RQIs, respectively. The rank n of an RQI is nothing 
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but the shortest period of its continued-fraction expansion, and the number of the members of the relevant MEC of 
RQIs is equal to n. 

The subject in the remaining part of this subsection concerns only in Sec. VII B in the text. Let 9 be an RQI 
specified by the ordered set (fco, fci, ■ ■ ■ , k n -t)- Then, the dual 9 to 9 is specified by the inversely ordered set. 
An MEC of RQIs is called self-dual or strongly self-dual if the corresponding MEC of normal UFMs is self-dual or 
strongly self-dual, respectively. A strongly self-dual MEC of RQIs includes an RQI 9 so that 9 = 9; the first period 
of the continued-fraction expansion of 9 is mirror-symmetric. We may say this RQI to be strongly self-dual. If £ is a 
member of a self-dual MEC of RQIs, £ is modular equivalent to £ or, equivalently, Z[£] and Z[£] are scaling-equivalent. 

A sufficient condition for an MEC of RQIs to be self-dual is that it includes an RQI £ so that £ + £ € Z or, 
equivalently, £ takes the form £ := p + q\f~D or £ := (p + qy/D ) /2, where p, q, and D are natural numbers. Therefore, 
the one period of the continued fraction expansion of a quadratic irrational of either of the two forms is mirror 
symmetric if it is considered to be a cycle. 



6. Quasi- reduced quadratic irrational 



The companion quadratic irrational £ of a quasi-normal UFM is a quasi-reduced quadratic irrational (QRQI) which 
is defined in Refs. [H^,[ll] because it satisfies the conditions: £ > 1 and £ > (or, equivalently, £ < 0). A necessary 
condition for a quadratic irrational being larger than unity to be an QRQI is that its continued-fraction expansion is 
periodic but for the first term. It can be shown readily that a QRQI £ is related to an RQI 9 as £ = 9 — k with k 
being an integer satisfying the inequality, < k < fco — 1 with fco := \ff\. 

The relation £ = 9 — k be twee n a QRQI £ and an RQI 9 is written in terms of the modular transformation specified 



by a UFM S given in Eqs. (All) as £ = S k (9). It follows from the equation 9 = M (9) with M being the companion 
matrix of 9 that £ is the companion qu adrat ic irrational of the quasi-normal UFM, Mk = S k MoS~ k : £ = Mfe(£). 
An explicit form of Mk is given by Eq. ( A16| ), where we have set Mq := (ko, ki, • • • , fcn-i)- The RQI 9 yields k, 



o 

QRQIs, £ = 9 — k with k = 0, 1, 2, • • • , fco — 1, which we shall call companion QRQIs of 9. Thus, there exists a 
bijection between the set of all the QRQIs and the set of all the quasi-normal UFMs. More strongly, there exists a 
bijection between the set of all the MECs of QRQIs and that of the quasi-normal UFMs. Both the two types of MECs 
combined by the bijection have structures of iV-cycles, and these structures are retained in the bijection. That is, the 
two iV-cycles are isomorphic. Therefore, we can specify an MEC of QRQIs by the symbol specifying the corresponding 
MEC of quasi-normal UFMs: ((fco, fci, ■ ■ ■ , fcn-i)) with n being the rank of the class. Thus, the relevant MEC of 
QRQIs includes n RQIs, 6j with j = 0, 1, 2, ■ ■ ■ , n—l and \9f\ = kj. The companion QRQIs of 9j are £^' > := 6j — fc 
with fc = 0, 1, 2, kj - 1. 

As a typical example of the MEC we take ((1, 2, 1)), which is composed of the four QRQIs: 

(0) _ l + VTo (d _ VTo 

50 — g > So — 1 < 2 ' 

(1) _ Vw (2) _ 2 + VTo 

51 - 2 , So - 3 ■ 
ie ratio r associated with this MEC is given b 

by 



The ratio r associated with this MEC is given by r = 3 + \/T0 • The companion matrix of £q° , for example, is given 



M =[ll)- (A28) 
This example can be used to check or understand various results appearing in the present paper. 
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FIG. 1. A derivation of a binary ID QL by the projection method. The mother lattice A is cut with the strip, and then 
projected onto Eu . 
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FIG. 2. A geometrical derivation of the IR of the Fibonacci lattice. The canonical set of basis vectors, {a, b}, is transformed 
by the hyperscaling operation T into another set, {a', b'}. The two types of new intervals, a 1 and are related to the original 
intervals, a and (3, by the IR. The window is shrunken by T from W to W' = t~ 1 W, while the sizes of the two types of intervals 
are expanded by the factor r. 
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FIG. 3. A geometrical illustration for the SR (3.2). This case is different from the one illustrated in Fig. ^ in that the third 
basis vector a + b is paired with a but not b to form a canonical set. 
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FIG. 4. The "finger print" of the MLD class ((1, 2, 1)). 
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